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I I Abstract: A successful implementation of thermal leptogenesis requires the re-heat tem- 

g_) ' perature after inflation Tr to exceed ~ 2 x 10^ GeV. Such a high Tr value typically leads 

, ■ to an overproduction of gravitinos in the early universe, which will cause conflicts, mainly 

with BBN constraints. Asaka and Yanagida (AY) have proposed that these two issues 
. can be reconciled in the context of the Peccei-Quinn augmented MSSM (PQMSSM) if one 

0^ , adopts a mass hierarchy m(sparticle) > m(gravitino) > m(axino), with m(axino) ~ keV. 

\ In this case, sparticle decays bypass the gravitino, and decay more quickly to the axino 

I LSP, thus avoiding the BBN constraints. In addition, thermally produced gravitinos decay 

CN ' inertly to axion+axino, also avoiding BBN constraints. We calculate the relic abundance of 

mixed axion/axino dark matter in the AY scenario, and investigate under what conditions 
a value of Tr sufficient for thermal leptogenesis can be generated. A high value of PQ 
breaking scale fa is needed to suppress overproduction of axinos, while a small vacuum 
^ , misalignment angle 9i is needed to suppress overproduction of axions. The large value of 

^ \ fa results in late decaying neutralinos. We show that, to avoid BBN constraints, the AY 

scenario requires a rather low thermal abundance of neutralinos, while higher values of 
neutralino mass also help. We combine these constraint calculations along with entropy 
production from late decaying saxions, and find the saxion needs to be typically at least 
several times heavier than the gravitino. A successful implementation of the AY scenario 
suggests that LHC should discover a spectrum of SUSY particles consistent with weak 
scale supergravity; that the apparent neutralino abundance is low; that a possible axion 
detection signal (probably with in the sub-//eV range) exists, but no direct or indirect 
signals for WIMP dark matter should be observed. 
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1. Introduction 

A wide assortment of data from atmospheric, solar, reactor and accelerator experiments 
can all be explained in terms of massive neutrinos with large mixing angles which undergo 
flavor oscillations upon propagation through matter or the vacuum [1]. The emerging 
picture of the physics behind neutrino oscillation data is most elegantly explained by the 
presence of massive gauge singlet right-hand Majorana neutrino states Ni [i = 1 — 3 a 
generation index) which give rise to see-saw neutrino masses [2]: m^.. ~ {f,y^v)^/MNi with 
fl,^ the neutrino Yukawa coupling, v the vev of the Higgs field, and M^. ~ 10^ — 10^^ GeV. 

In addition to explaining neutrino oscillation data, the presence of massive Ni states of- 
fers an elegant explanation of baryogenesis in terms of leptogenesis [3] , wherein the states 
Ni exist in thermal equilibrium in the early universe, but decay asymmetrically to lep- 
tons versus anti-leptons. The lepton-anti-lepton asymmetry is then converted to a baryon 
asymmetry via B and L violating, but B — L conserving, sphaleron effects [4] . To realize 
the thermal leptogenesis scenario, the lightest of the heavy neutrino masses Mi must be 
> 2 X 10^ GeV. In order to produce such states thermally, a re-heat temperature of the 
universe after inflation of Tr > Mi > 2 x 10^ GeV is required [5] . 

Augmenting the Standard Model with a new, extremely high energy scale Mf^^ natu- 
rally leads to severe quadratic divergences in the Higgs sector which will need to be highly 
fine-tuned. The solution here is to also incorporate supersymmetry (SUSY), which reduces 
quadratic divergences to merely logarithmic, and ameliorates the fine-tuning problem [6]. 
While the addition of weak scale softly broken SUSY into the theory is actually supported 
by the measured values of the gauge couplings from LEP experiments, it also leads to 
new conundrums such as the gravitino problem: the production of gravitinos in the early 
universe can lead to (i) overproduction of LSP dark matter {e.g. the lightest neutralino) 
beyond relic density limits obtained from WMAP and other experiments, or (ii) disrup- 
tion of the successful explanation of Big Bang nucleosynthesis by introducing late decaying 
quasi-stable particles whose decay products can break up the newly minted light elements. 
The common solution to the gravitino problem [7] is to require a sufficiently low re-heat 
temperature such that thermal gravitino production is suppressed enough to avoid over- 
production of dark matter or disruption of BBN [8] . For gravitino masses in the few TeV 
or below range, a value of Tr < 10^ GeV is required. Naively, this is in obvious conflict 
with the Tr requirements of thermal leptogenesis. 

A variety of solutions have been proposed to reconcile leptogenesis with the SUSY 
gravitino problem. One is to abandon the "thermal" aspect of leptogenesis, and invoke 
non-thermal leptogenesis wherein the heavy neutrino states are produced via some scalar 
field decay, for instance the inflaton [9]. Another suggestion is to invoke the gravitino as 
LSP, so it does not decay. However, the gravitino LSP scenarios fall back into the BBN 
problem since then the NLSP SUSY particle suffers a late decay into gravitino plus SM 
states which again injects high energy particles into the post-BBN plasma. One solution 
is to speed up NLSP decay via a small component of i?-parity violation [10, 11]. 

In a recent work [12], we proposed an alternative scenario, invoking mixed axion/axino 
dark matter, as would occur in the Peccei-Quinn [13-16] augmented MSSM (the PQMSSM) 
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[17,18]. In this case, we invoked models with very heavy gravitinos, mg > 10 TeV, so that 
gravitinos decay before the onset of BBN. Then, overproduction of dark matter can be 
avoided by requiring an axino LSP with mass rria ~ 0.1 — 1 MeV. Neutralinos produced 
either thermally or via gravitino decay will themselves decay typically to states such as dj, 
so that the dark matter abundance is reduced by a factor rria/m^^ [26]. The bulk of dark 
matter then resides in thermally produced axinos and/or in axions produced from vacuum 
mis-alignment. By driving up the value of PQ breaking scale fa/N, thermal production of 
axinos is suppressed, and higher values of Tr are required to maintain a total axino plus 
axion relic abundance of f^aa^^ ~ 0.11. To avoid overproduction of axions at high fa/^, 
we adopted a small vacuum mis-alignment angle 6i ~ 0.05. However, the large values of 
fa/N ~ 10^^ — 10^'^ GeV suppress the Zi decay rate, thus interfering with BBN from a 
different avenue. Models with a high-mass, bino-like Zi and low "apparent" ^'^-^h? can 
avoid the BBN bounds, and allow Tr values in excess of 10^*^ GeV to be attained. As we 
showed, such conditions with ~ 10 — 30 TeV can be reached in Effective SUSY [19, 20] 
or mirage-unification SUSY breaking [21] models. 

A related scenario to reconcile thermal leptogenesis with the gravitino problem — 
using mixed axion/axino dark matter — was proposed much earlier by Asaka and Yanagida 
(AY) [22]. Their proposal was to work within the context of the PQMSSM, but with a 
sparticle mass hierarchy m(sparticle) > > m^. In this case, the couplings of MSSM 
sparticles to axinos are larger than the couplings to gravitinos, so that the long-lived decays 
to gravitino are bypassed, and the sparticles are assumed to decay to an axino LSP shortly 
before the onset of BBN. Furthermore, thermally produced gravitinos decay inertly via 
G — 7- aa and so do not disrupt BBN. Reheat temperatures as high as Tr ~ 10^^ were 
claimed to be possible. 

In this paper, we re-visit the AY scenario, incorporating several improvements into our 
analysis. In particular, we implement 

1. the latest astrophysically measured value of [23] 

Odm/i^ = 0.1123 ± 0.0035 at 68% CL; (1.1) 

2. the latest calculations for thermal production of gravitinos and axinos; 

3. vacuum-misalignment production of axions as an element of the dark matter abun- 
dance; 

4. the latest BBN constraints on late decaying particles; and finally 

5. a careful treatment of entropy production from late decaying saxions. Since entropy 
production from saxion decay will also dilute the matter-antimatter asymmetry by a 
factor r (to be defined later), in this case a re- heat temperature Tr > 2r x 10^ GeV 
will be needed. 

The re-analysis of the AY scenario taking into account points 1.-4. is presented in Sec. |2|, 
while entropy injection from saxion decay is discussed in detail in Sec. ^. In Sec. |^, we 
present our final conclusions and consequences of the AY scenario for LHC physics and 
dark matter searches. 
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2. Relic density of mixed axion/axino DM in the AY scenario 



2.1 MSSM parameters 

The only relevant MSSM parameters for our analysis are the Z\ and G masses and 
mg, the Z\ bino component v'^ in the notation of [6], and the Z\ abundance after freeze- 
out ^2,^- '^'^^ remaining of the MSSM parameters only impact the running of the gauge 
couplings and the value of fig, which depend on all the gaugino masses (see Eq. ( |2.5D below). 
However, as shown below, in the AY scenario with Tr > 10^ GeV, the contribution from 
G — )• aa decay to the dark matter relic density is negligible. Thus the dependence on the 
entire SUSY spectrum is very mild. 

For illustration we will use a generic mSUGRA scenario with mg = 1000 GeV, = 
1000 GeV, Ao = 0, tan/3 = 55 and /i > 0, which gives m^^= 430 GeV and f^^i = ^-^^^ 
treat and as free parameters throughout most of the numerical analysis. The bino 
component of the Z\ wavefunction is important, since it determines the Z\ — a coupling. 
For simplicity we will assume a purely bino neutralino, which is valid for a large portion 
of the mSUGRA parameter space. We also take = ra-^jl and ma < mg in order to 
have an axino LSP with a gravitino NLSP. 

2.2 Mixed axion/axino abundance calculation 

Here, we consider four mechanisms for dark matter production in the AY scenario. 
• If the reheat temperature Tr exceeds the axino decoupling temperature 



axinos will be in thermal equilibrium, with an abundance given by 

^l^'h' ^ 0.38 { . (2.2) 



1 keV 



To avoid overproducing axino dark matter, the RTW bound [25] then implies that 
m-a < 0.3 keV. 

In the case where Tr < T^cp, the axinos are never in thermal equilibrium in the early 
universe. However, they can still be produced thermally via radiation off of particles 
that are in thermal equilibrium [27,28]. Here, we adopt a recent calculation of the 
thermally produced (TP) axino abundance from Strumia [29]: 

nj^h^ = 1.24gtF(g3) ( (2.3) 

ya ^y'^^GeVlO^ GeV V/a/iVy 

with -F(g3) ~ 20(7| In and is the strong coupling constant evaluated ai Q = Tr. 

In super symmetric scenarios with a quasi-stable neutralino, the Zis will be present 
in thermal equilibrium in the early universe, and will freeze out when the expansion 
rate exceeds their interaction rate, at a temperature roughly Tjo ~ The 
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present day abundance can be evaluated by integrating the Boltzmann equation. 
Several computer codes are available for this computation. Here we use the code 
IsaReD [30], a part of the Isajet/Isatools package [31,32]. 

In our case, each neutralino will undergo decay to the stable axino LSP, via decays 
such as Z\ — >■ 07. Neutralinos may also decay via e.g. Z\ — >■ G7, but these modes 
are suppressed by rather than l/{falNY , and so the decay to gravitinos is 

suppressed (one of the hallmarks of the AY scenario). Thus, the non-thermally pro- 
duced (NTP) axinos will inherit the thermally produced neutralino number density, 
and we will simply have [27] 

• Since here we are attempting to generate reheat temperatures Tr > 10^ GeV, we 
must also include in our calculations the thermal production of gravitinos in the 
early universe. Wc adopt the calculation of Pradler and Steffen in Ref. [33], who 
have estimated the thermal gravitino production abundance as 

(2.5) 

where Ui = (0.018,0.044,0.117), h = (1.266, 1.312, 1.271), gi are the gauge couplings 
evaluated at Q = Tr and Mj are the gaugino masses also evaluated at Q = Tr. For 
the temperatures we are interested in, this agrees within a factor of about 2 with 
the more recent calculation by Rychkov and Strumia [33], which is sufficient for our 
purposes. 

Since the only kinematically allowed gravitino decay mode is to an axion plus an 
axino LSP, the abundance of axinos from gravitino production is given by 

nih' = ^m'h\ (2.6) 

rriQ ^ 

while the abundance of axions from gravitino production is given by 

nGh'' = ^n^/h". (2.7) 

niQ G 

For axino masses in the MeV range and gravitino masses in the TeV range, the 
prefactor above is extremely small, making the contribution from gravitino decays to 
the axino relic abundance negligible, what allows us to evade overproduction of dark 
matter via thermal gravitino production. 

• Here, we consider the scenario where the PQ symmetry breaks before the end of 
inflation, so that a nearly uniform value of the axion field 6i = a{x)/{fa/N) is ex- 
pected throughout the universe. The axion field equation of motion implies that 
the axion field stays relatively constant until temperatures approach the QCD scale 
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Tqcd ~ 1 GeV. At this point, the temperature-dependent axion mass term turns on, 
and a potential is induced for the axion field. The axion field rolls towards its mini- 
mum and oscillates, filling the universe with low energy (cold) axions. The expected 
axion relic density via this vacuum mis- alignment mechanism is given by [34] 

n^h'c.0.23m)e^^j^l^y (2.8) 

where < Oi < n and f{0i) is the so-called anharmonicity factor. Visinelli and 
Gondolo [34] parametrize the latter as f{Oi) = In yjz^J^J ■ '^^^ uncertainty in 
Qah'^ from vacuum mis-alignment is estimated as plus-or-minus a factor of three. 

In this paper, we will evaluate the mixed axion/axino relic density from the above four 
sources: 

Qaah"^ = ^Fh^ + ^ih^ + ^^h^ + ^^h^ + fiah^- (2.9) 

Over much of parameter space, if nia is taken to be of order the MeV scale or below, then 
the contributions $7?', Q,?^ and are subdominant. 

In Fig. I, we illustrate in frame a) the relative importance of the four individual 
contributions as a function of fa/N, for ^^_^h'^ = 10, m^^ = 430 GeV and = 0.5m^^. 
For the axion/axino sector we take 9i = 0.05 and rria = 100 keV. The value of Tr is adjusted 
such that ^aah^ = 0.1123. Low values of 6i suppress axion production and allow higher 
values of fa/N to be probed; the higher values of /a/-/V suppress thermal axino production, 
thus allowing for higher Tr values to compensate. 

For low fa/N values, the TP axino contribution is dominant. But as fa/N increases, 
the axion component grows. For higher values of fa/N, the vacuum-misalignment produced 
axion component dominates, and the dark matter is predominantly composed of cold ax- 
ions. The contribution of axino dark matter from Zi and G decays are always negligible 
in this case. 

In frame b) of Fig. |l], we show the value of Tr which is needed to enforce the total abun- 
dance of mixed axion/axino dark matter to be ^aah"^ = 0.1123. We show cases for = 0.1 
and 1 MeV. As /a/-/V increases, the axino-matter coupling decreases, and one would expect 
the thermal production of axinos to decrease. Since we enforce O.aah'^ = 0.1123, then higher 
values of Tr are needed to compensate and enhance the thermal production of axinos [35] 
(and gravitinos). We see that the value of Tr can be pushed to over 10^ GeV for rria = 1 
MeV, and to over 10^*^ GeV for rria = 0.1 MeV, thus allowing high enough Tr for thermal 
leptogenesis. 

2.3 Constraints from cold/warm dark matter 

Depending on its mass, the axino might constitute warm (WDM) or hot (HDM) dark 
matter; the latter possibilities are severely constrained by the matter power spectrum and 
reionization [27,36] (see also [37,38]). We consider axinos with mass 1 — 100 keV as mostly 
WDM, and axinos with mass < 1 keV as mostly HDM. Since these bounds on the amount 
of WDM/HDM are model dependent [36], we do not impose strict WDM/HDM constraints 
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Figure 1: Upper frame: Contribution of axions and TP and NTP axinos to the DM density as 
a function of the PQ breaking scale fa/N, for m^^ = 430 GeV, ^^i^^ ~ ™s — 100 keV and 
9i = 0.05; Tr is adjusted such that ftad.h^ — 0.1123. We assume = Lower frame: the 

value of Tr that is needed to achieve ilaah^ = 0.1123 for ma = 0.1 and 1 MeV. 



on our results. However, for guidance, we will keep track of PQMSSM parameter points 
with potentially too large WDM and/or HDM components^: As in [12], we disfavour points 
with 

^a/^aa > 0.2 VI keV <ma< 100 keV, (WDM) 
Qa/^aa > 0.01 V < 1 keV, (HDM) 



(2.10) 



where 0^ = 0,7^ + 0^ + Of. This is rather conservative. A rough estimate based on the 
neutrino mass limit [38] from cosmological data, ^ rrii, < 0.41 to 0.44 eV, gives that up 
to 4-5% HDM contribution could be acceptable. Moreover, Boyarsky et al. in [37] found 
that in case of a thermal relic (TR), 100% WDM is allowed for mxR > 1.7 keV, while for 



^Axions produced from gravitino decay will also constitute HDM. However, since this contribution to 
the total DM density is suppressed by rria/mg, it can be safely neglected. 
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iT^TR = 1-1 keV as much as 40% WDM is allowed at 95% CL. We will also indicate these 
bounds, which are considerably weaker than Eq. ( 2.10| ). 



2.4 Constraints on Zi decay from BBN 

The AY scenario naturally avoids BBN constraints on late decaying gravitinos by assuming 
the mass relation m(sparticle) > > ma, so that the G decays inertly 100% of the time 
into aa. However, by searching for PQMSSM parameters which allow Tr ^ 2 x 10^ GeV 
while avoiding overproduction of mixed axino/axion dark matter (the latter requires large 
fa/N ~ 10^^ GeV and small 9i), we have pushed the Zi lifetime uncomfortably high, so 
that its hadronic decays in the early universe now have the potential to disrupt BBN. The 
Zi lifetime and hadronic branching fraction is calculated in Ref. [12,27]. 

Constraints from BBN on hadronic decays of long-lived neutral particles in the early 
universe have been calculated in Ref 's [39-41]. Here, we adopt the results from Jedamzik [41] 
The BBN constraints arise due to injection of high energy hadronic particles into the ther- 
mal plasma during or after BBN. The constraints depend on three main factors: 

• The abundance of the long-lived neutral particles. In Ref. [41], this is given by VLxh?' 
where X is the long-lived neutral particle which undergoes hadronic decays. In our 
case, where the long-lived particle is the lightest neutralino which decays to an axino 
LSP, this is just given by the usual thermal neutralino abundance /i^, as calculated 
by IsaReD [30]. 

• The lifetime tx of the long-lived neutral particle. The longer-lived X is, the greater 
its potential to disrupt the successful BBN calculations. In our case tx = r^^ oc 
{fjNf/m\^. 

• The hadronic branching fraction of the long-lived neutral particle. If this is very 
small, then very little hadronic energy will be injected, and hence the constraints 
should be more mild. 

From the above list, we see that BBN directly constrains the MSSM parameters and 
m^^) of the PQMSSM model. The constraints also depend on /q, since its value directly 
affects T^^. The BBN constraints are shown in Fig. 9 (for mx = 1 TeV) and Fig. 10 (for 
mx = 100 GeV) of Ref. [41], as contours in the tx vs. Qxh'^ plane, with numerous contours 
for differing values ranging from 10^^ to 1. For Bf^ ~ 0.1, for instance, and very large 
values of Qxh'^ ~ 10 — 10^, the lifetime tx must be < 0.1 sec, or else the primordial 
abundance of ^He is disrupted. If ilxh? drops below ~ 1, then much larger values of tx up 
to ~ 100 sec are allowed. If one desires a long-lived hadronically decaying particle in the 
early universe with tx ^ 100 sec, then typically much lower values of ^Ixh? ~ 10^^ — 10~^ 
are required. 

We have digitized the constraints of Ref. [41], implementing extrapolations for cases 
intermediate between values of parameters shown, so as to approximately apply the BBN 
constraints to the AY scenario with a long-lived neutralino decaying during BBN. The 
results are shown in Fig. ^, as a function oim^^. Prom Fig. ^b), we have Tr > 10^ GeV for 
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Figure 2: BBN bounds on late decaying neutralinos {Zi — > Z/j + a) from Ref. [41] as a function of 

for different values of fa/N , assuming a bino Zi and 3> r 
the curves are excluded by BBN constraints. 



The values of ft?; above 
^1 



fa/N ~ 10^3 GeV. From Fig. |, we see that, for m^^ ~ 400 GeV and fa/N = 10^^ GeV, 
we need Vi^ h? < 0.4 in order to satisfy the BBN bounds. In particular, the assumed value 
for ri^^/i^ (=10) in Fig. | only satisfies the BBN constraints for fa/N < 2 x 10^^ QgV or 
^-R ^ 3 X 10*^ GeV. Therefore we see that the AY scenario with Tr > 10^ GeV requires 
quite low values of ^^_^h?' and/or a heavy Zi. 

2.5 Scan over PQ parameters 

The results of the last section were for a specific choice of the axino mass and 6i values. 
Next, we examine which values of are possible for arbitrary values of /a, rria and di. For 
now we keep the MSSM parameters {Q^^Ji^ and wi^^) fixed. As discussed in the previous 
section, the BBN constraints in general require low ^^Jt^ and high ni^^. We therefore 
assume = 0.04 and = 430 GeV, which are values consistent with, e.g., an 

mSUGRA point near the apex of the Higgs funnel region. To probe the full PQ parameter 
space we perform a random scan over the PQ parameters in the range 



ma G [10~^ 10] GeV, 
fa/N G [10^, 10^^] GeV, 
e^ e [0, tt] . 



(2.11) 



and calculate the value of Tr which is needed to enforce ^aah^ = 0.1123. The results of 
our scan are shown in Fig. |3|, where we plot the derived value of Tr versus PQ breaking 
scale fa/N. In the plot, dark blue and dark red points have mainly CDM with at most 



20% WDM and at most 1% HDM admixture, c.f. Section |j. The hght blue and light red 
points have higher values of WDM or HDM. The red points are excluded by bounds derived 
from Ref. [41] on late decaying neutralinos which could destroy the succesful predictions 
of BBN. Blue points are allowed by BBN constraints. Applying the WDM bounds of 
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nig =430GeV,Qg h^ = 0.04 




fJN(GeV) 

Figure 3: Allowed and disallowed points in the fa vs. plane for m^^ = 430 GeV and ^^_h^ = 
0.04, including BBN constraints on late Zi decay. Dark blue points have mainly CDM with at most 
20% WDM and at most 1% HDM admixture; the dashed blue line indicates the WDM limit by 
Boyarsky et al. [37] and up to 5% HDM 



Boyarsky et al. [37] for a thermal relic and allowing up to 5% HDM, the boundary between 
dark and light blue points moves left to the dashed blue line. 

We do see that a number of dark blue points with mainly CDM, and which also respect 
BBN bounds, can be generated with Tr > 2 x 10^ GeV. These thermal leptogenesis- 
consistent points all require fa/N > 10^^ GeV (or fa/N > 6 x 10^^ GeV with weaker 
WDM/HDM requirements). 

In Fig. we show the same scan in the 9i vs. Tr plane. Here, we see that the 
CDM/BBN consistent points with high Tr all need rather small values of axion mis- 
alignment angle 6i < 0.5. This is needed since, at large Tr, large fa/N is necessary to 
suppress overproduction of axinos, while small 6i is needed to suppress over-production of 
axions. 

In Fig. |5|, we show the same PQMSSM parameter scan for m^^ = 430 GeV and ^^J^"^ = 
0.04, but in the Q-ah^ vs. Tr plane. Here, we see that the CDM/BBN consistent points 
with high Tr can have both large and small values of Q.ah? . Solutions with Oq/i^ ~ 0.1 
usually have very light axinos (to suppress Oq/i^) and moderate 9i values. Solutions with 
< 0.1 usually have heavier axinos and small 9i. 

While the results of the previous figures are restricted to specific values of and 
Cl^^h"^, the overall scheme is much more general. The lesson here is that the AY scenario 
for reconciling thermal leptogenesis with the gravitino problem can work provided certain 
conditions on SUSY models are met. These conditions are rather similar to those needed 
by the large scenario put forth in Ref. [12]. After adopting a model with a sparticle 
mass hierarchy of m(sparticles) > > rria, with ma ~MeV scale and Mweak, one 

needs the following features: 

• To allow for Tr > 2 x 10^ GeV, one must suppress thermal production of axinos via 
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Figure 4: Allowed and disallowed points in the 9i vs. Tr plane for = 430 GeV and f^^i^^ 
0.04, including BEN constraints on late Zi decay. 
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Figure 5: Allowed and disallowed points in the Q.ah?' vs. Tr plane for m^^ = 430 GeV and 
ris — 0.04, including BBN constraints on late Zi decay. 



a large value of fa/N > lO^^ Q^y 

• To suppress overproduction of axions one must adopt a lower range of mis-alignment 
angle 9i < 0.5 (or 6i < 0.8 taking into account the factor 3 uncertainty in Eq. ( |2.8D ) . 

• The large value of fa/N increases the Zi lifetime, which brings in BBN constraints 
on late-decaying neutral particles. To avoid BBN bounds, it helps to invoke 1. a bino- 
like Zi so that 1, 2. a low apparent neutralino relic abundance ^'^-^h? < 1 and 

3. a large value of to help suppress the Zi lifetime. 
These conditions are illustrated in a more model independent way in Fig. |6|. Here, we 
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Figure 6: Allowed and disallowed points in the ^^i^^ vs.m^^ plane for a general scan over SUSY 
models with a bino Zi. For all points, we require Tr > 2 x 10^ GeV and assume = m^^/2. 
Dark blue points are consistent with BBN and have mainly CDM with at most 20% WDM and/or 
1% HDM admixture. The region below the dashed line represents the MSSM parameter space where 
99% of the DM/BBN consistent solutions lie when applying weaker WDM/HDM requirements as 
discussed in the text. 

assume the AY mass hierarchy, but extend our previous scan to the whole PQMSSM 
parameter space: 

m-a e [10^^ 10] GeV, 
fa/N G [10^, 10^5] GeV, 
Oi G [0, ^] , 

n^y e [lO"^lO^], 

m^^ e [10, 10^] GeV . 

As before, we assume = m^j2 and the blue points are BBN-allowed, while red points 
violate BBN bounds. The dashed line indicates the boundary blow which 99% of the 
DM/BBN consistent solutions lie when applying weaker WDM/HDM requirements (WDM 
limit according to Boyarsky et al. [37] and up to 5% HDM, cf. Fig. |3|) This line can be 
interpreted as a natural upper bound for as a function of "m^^. From this we see 

that models with m^^ < 10 GeV require ^ 10^^, while values of as high as 

10'^ can be consistent with thermal leptogenesis if the neutralino is in the TeV range. 

3. Dilution of DM by entropy production from saxion decay 

Up to this point, we have neglected an important element of the axion supermultiplet, 
namely the spin-0 saxion field s{x) which is expected to obtain a soft SUSY breaking mass 
at the GUT scale of nis ^ rriQ. In the same way as axinos, saxions can be produced in 
thermal equilibrium (if Tr > T^cp) or out of equilibrium from scatterings of particles in 
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the plasma (if Tr < Td^p). However a second mechanism of saxion production is also 
possible [45]. After supersymmetry breaking, the saxion potential V{s) develops a global 
minimum, causing the saxion field to coherently oscillate around its minimum. This coher- 
ent oscillation can have a large energy density, which contributes to the total saxion energy 
density if Tr < T^cp- However, if Tr > T^cp, the coherent oscillating saxions will couple to 
the thermal plasma and thermalize. 

Once the saxions decouple from the thermal plasma (at T = T^cp) and become non- 
relativistic (at T ~ rus), their energy density (ps) scales as (or R~^), while the thermal 
plasma's energy density (prad) scales as (or If the saxion lifetime is sufficiently 

long, at some temperature Tg, we will have ps > Prad a-i^d the universe will become tem- 
porarily matter dominated until the saxions decay. 

Being a i?-parity even state, the saxion can decay to standard model states or pairs 
of sparticles. Since we assume ~ mo, the decay into SUSY states will be kinematically 
suppressed and the saxion decays will mostly consist of SM particles. Therefore the saxion 
decay products will thermalize in the thermal plasma, which is then "relatively reheated" 
[47] with respect to other decoupled particles, such as axinos. As a consequence, all particles 
decoupled from the thermal plasma during the saxion decay will have their number density 
diluted with respect to the thermal bath's. Below we introduce the relevant expressions 
necessary for computing this dilution factor (r) and in Sec. we discuss how the inclusion 
of the saxion field impacts our previous results. Since entropy injection from saxion decay 
will also dilute the matter-antimatter asymmetry by a factor r, in this case a re- heat 
temperature Tr > 2r x 10^ GeV will be required for a successful implementation of thermal 
leptogenesis. 

3.1 Saxion production and decay 

As mentioned above, if Tr exceeds T^cp in the early universe, saxions are produced in 
thermal equilibrium such that 

y.m. = ^ 10-3^, (3.1) 
s GeV 

where s = 27r^(7*T^/45 is the plasma entropy density and is the saxion yield. For T/j < 
Tdcpi saxions can still be produced thermally, although to our knowledge a full calculation 
is not yet available. In Ref's [22] and [45], thermal saxion production is estimated to be 

which we will adopt for our calculations. 

In addition, saxions can be produced via coherent oscillations of the saxion field in the 
early universe. Although the exact mechanism depends on the saxion potential near the 
SUSY breaking scale, the energy density associated with the coherent oscillations can be 
parametrized by the initial saxion field strength, Sj. Natural values for Sj are ~ fa or 
Si ~ Mpi. Unless stated otherwise, we will assume Si = fa/N. The saxion energy density 
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Figure 7: Saxion yield Ys versus fa/N for Tji — 10^ GeV and rris = 0.1 and 1 TeV. We assume 

S^/{fa/N) = 1. 



is estimated for the case of very high values of Tr with Tj > (here, Tr is related to the 
inflaton decay width Tj as Tr = (S/g^vr^) V4(j\^^^r^)i/2) [45] 

m., \ 1/2 ^ ""^^ ^ ^ 



— ~ 1.5 X 10" 

s 



GeV 

while for Tj < rrig, 

— ~ 2.1 X 10-9 GeV 



1 GeV; 



105 GeV / 



ifa/N) Y 

1012 GeV J 



Ua/N) 



ifa/N) Y 

1012 GeV J 



(3.3) 



(3.4) 



The summed saxion abundance is then given by the thermal production if Tr > T^cp or by 
the sum of thermal production plus the abundance from coherent oscillations, if Tr < T^cp- 

As an example, we show in Fig. ^ the saxion yield Yg versus fa/N for = 0.1 and 
1 TeV and for Tr = 10^ GeV. At low fa/N, Tr > Tdcp, and saxions are produced in 
thermal equilibrium. Once T^cp rises above Tr, thermal saxion production dominates, but 
decreases with increasing fa/N until the point where saxion production from coherent 
oscilations dominates. 

The saxion is an i?-parity even state which is expected to dominantly decay into two 
gluons: s gg. The saxion decay width differs by factors of two in Ref's [46], [22] and [45]. 
By an independent calculation, we find 



r(s ^ gg) 



1 s 



(3.5) 



327r3(/JiV)2^ 

in agreement with [22]. The saxion may also decay (or not, model dependently) via s — >■ aa, 
and in the DFSZ [16] model, into qq or ££. These latter decays are suppressed in the KSVZ 
model [15]. Saxion may also decay to ZiZj, 77 and gg. For saxion decay to gluino pairs, 
we calculate the interaction as2 

as 



C 3 



Mfa/N) 



(3.6) 



^Our saxion-gluino-gluino interaction differs by a factor of 2 from Ref. [29]. 
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Figure 8: Decay widths for s ^ gg and s ^ gg a,s a function of nig for fa/N — 10^^ GeV and 
TO, = 1 TeV. 



and then find 



r(s ^ gg) 



4m 



2\ 3/2 



(3.7) 



The two widths are compared in Fig. |8|, where the s ^ gg decay is found to always 
dominate. 

The temperature associated with saxion decay and entropy injection is given by [47] 



(3.8) 



All the saxion decays essentially finish entropy injection by the time the universe cools to 
this value [47]. For simplicity, here we will assume the F^ = F(s — ?> gg) so that our results 
are independent of nig. Folding in the additional strong decay s ^ gg will result typically 
in a small increase in T^. In Fig. ^, we plot the value Tg as a function of fa/N for three 
different values of the saxion mass. 

3.2 Entropy injection from saxion decay 

Armed with expressions for the saxion production rate in the early universe, we next 
calculate the temperature Tg at which the saxion energy density in the universe equals the 
plasma energy density: 

(3.9) 



Ps{Te) = Prad{Te) = ^^-^T^. 



30 



Using ps = mgYgS and s = '^^g^T^ , we find 



(3.10) 



If Tg exceeds (i.e. if saxion domination occurs before saxion decay), then saxions 
dominate the energy density of the universe for < T < Tg. In this case, saxion decays 
may inject significant entropy and dilute whatever abundances are present at temperature 
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Ts . The situation is shown in Fig. |To|, where we show the value of (blue horizontal lines) 
and the value of Tg (red lines) for fa/N = 10^^ and 10^^ GeV. For Tr greater than the 
r = 1 intersection points, significant entropy injection can occur. 

If r > 1, Ts must satisfy Tg > 10 MeV, so the universe becomes radiation dominated 
before BBN starts. On the other hand, if r < 1, the early universe is always radiation 
dominated and the usual BBN constraints on late decaying particles can be applied [41], 
as in the case of neutralino decays. If the conditions 10 MeV < Tg and Tg > Tg hold, 
then saxion decay can inject substantial entropy and dilute whatever relics are present and 
decoupled from the thermal plasma at the time of saxion decay (T = Tg). The ratio of 
entropy injection before and after a quasi-stable particle decay, for a matter dominated 
universe (Tg > T^), has been calculated in Scherrer and Turner [47], and is given by 

Sf 1/4 Ysfng Tp 
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where 5* is the number of relativistic degrees of freedom averaged over the saxion decay 
period, which we approximate by g^:(Ts). The above expression for r is only vahd for 
Tg > Tg (saxion dominated universe) or r > 1. However, if the saxion energy density never 
dominates the universe, the entropy injection is neghgible [47]. Therefore we assume r = 1 
(no entropy injection), if Tg < Tg. 

Assuming Tg is dominated by the s ^ gg decay, we plot in Fig. |ll|a) the value of r 
in the (fa/N) vs. Tji plane for rus = 0.1, 1 and 10 TeV, assuming Si/{fa/N) = 1 (for 
production from coherent oscillations). The solid lines all maintain Tg > 10 MeV, while 
dashed lines violate this constraint. We see first that for rrig = 100 GeV, significant entropy 
production only occurs for fa/N < 4 x 10^^ GeV; for larger fa/N, Tg is suppressed and the 
saxion lives long enough to decay during or after BBN. For rrig = 1 TeV, entropy injection 
can occur for fa/N < 10^^ GeV. 

If T/j lies below the r = 1 contours, then not much entropy is injected, but for high 
Tr, large entropy injection is possible and must be accounted for. The various contours of 
constant r initially increase with Tr. In this case, the saxion production is dominantly ther- 
mal. When the curves turn over, saxion production is dominated by coherent oscillations. 
In this case, as fa/N increases, the saxion field strength also increases (since Si/{fa/N) 
is fixed to 1), and much lower Tr values are allowed for substantial entropy production. 
Another noteworthy feature is that the contours of entropy production increase with Tr 
as nis increases. Thus, the dilution of DM from saxion decay can be reduced by requiring 
rather heavy saxions. Finally, when we compare Fig. ^ to Fig. 11, we see that the range 
of Tr ~ 10^ — 10^^ GeV for fa/N ~ 10^^ — 10^'^ implies that entropy dilution from saxion 
decay needs to be accounted for in our calculations for the case where Si/{fa/N) ~ 1 and 
nig = niQ = 1 TeV. 



In Fig. 11&), we plot again the entropy ratio contours, but this time taking Si/{fa/N) = 
0.1. In this case, saxion production from coherent oscillations is suppressed by the smaller 
initial saxion field strength value. This expands the range of large Tr at high fa/N where 
entropy injection is negligible. In cases such as these, the results of the previous section 
(and also Ref. [12]) remain viable, and entropy injection from saxion decay would be a 
negligible effect. From here on, we will assume Si/{fa/N) = 1. 

As mentioned before, the entropy injection from late decaying saxions will dilute the 
number density of any particle decoupled from the thermal plasma at T = Tg. Therefore, 
depending on Tg, the saxion production and decay may dilute thermally produced axinos, 
gravitinos, the quasi-stable Zis, sometimes the axions and of course the matter-antimatter 
asymmetry itself. In this latter case, the baryon-to-photon ratio r/^ = ns/n^ is diluted to 
a value ris/r. 

To accommodate the latter possibility, we note that the lepton asymmetry and con- 
sequently the baryon assymetry is proportional to the mass of the lightest right-handed 
neutrino. Mi (assuming Mi ^ M2^z) [48]. Therefore, to compensate the saxion dilution, 
heavier neutrinos are necessary. Since thermal leptogenesis requires Tr > Mi, in the case 
of saxion entropy injection, we need Tr, > r x (2 x 10^) GeV. 

To include the above effects into our new analysis, we adopt the following procedure: 
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Figure 11: Ratio of entropy r before and after saxion decay in the fa/N vs. Tji plane for mg = 
0.1, 1, 10 TeV and for a) Si/{fa/N) = 1 and h) Si/{fa/N) = 0.1. The dashed hnes correspond 
to Tdecay < 10 McV, whcn the entropy from saxion decay is injected after the beginning of BBN; 
these regions are Ukely excluded. 



• Calculate the thermal plus coherent oscillation yield of saxions 1^ in the early uni- 
verse. 

• Calculate the saxion decay temperature Tg. 

• Calculate Te and determine if saxions can dominate the universe (Te > Tg). 

• Calculate the final/initial entropy ratio r. 

• If r < 1, then the entropy injection is negligible and we require the saxion lifetime 
and relic density to satisfy the BBN bounds for late decaying particles in a radiation 
dominated universe, 
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• liTg < 10 MeV and r > 1, then the point is excluded due to entropy injection during 
or after BBN. 



• If T, > 10 MeV and r > 1, 

— dilute thermally produced axinos by factor r. 

— dilute thermally produced gravitinos by r. 

— If Ts < Tqcd = 1 GeV, dilute mis- alignment produced axions by r. 

— If Tg < Tfo = m^y25, dilute quasi-stable neutralinos by r. This condition can 
dilute axinos produced by neutralino decay, but also impacts the quasi-stable 
neutralino BBN bounds from Fig. 0. 



Our first results are shown in Fig. 12 for the same PQMSSM parameters used in 
Fig. but including saxion entropy injection with iris = 1 TeV. In frame a), we plot 
the relic abundance of thermally produced axinos (red), axions (blue), gravitino produced 
axinos (lavender) and neutralino produced axinos (magenta). The value of Tr is always 
adjusted to maintain Qaa = 0.1123, and so T/j/r is shown in frame b) for = 0.1 and 1 
MeV. 

For low values oi fa/N, the relic abundance curves track the values shown in Fig. ^. In 
this case, Tr is much lower than the leptogenesis value 2 x 10^ GeV, and the thermal yield 
of saxions is too low for significant entropy production. As fa/N increases, the thermal 
axino production drops, and the value of Tr must compensate by increasing the thermal 
yield of axinos so that ^adh^ = 0.1123 is maintained. At around fa/N ~ 10^^ GeV, the 
value of Ts drops below Tg, and significant entropy production from saxion decay occurs. 
The entropy injection dilutes the thermal axino and also axion production, so that a sharp 
increase in Tr is needed to offset the dilution effect: the dark matter abundance remains 
dominated by thermal axino production. However, the axion abundance is independent of 
Tr, and so its dilution due to saxion decay is plain to see in frame a). Since the entropy 
injection from saxion decay also dilutes the matter-antimatter asymmetry, we also show 
the trajectory of Tji/r once entropy injection is started. 

The value of Tr needed to maintain ^aah^ = 0.1123 increases sharply until the regime 
Tr > fa/N is reached. For such high values of Tr, the PQ symmetry is restored during 
re-heat, and re-broken during subsequent cooling. The universe should break into domains 
of different 9i and Si values (see e.g. M. Turner in Ref. [34]), and a modified treatment of 
dark matter will be needed. Therefore we neglect such solutions and impose the condition 
Tr < fa/N to our solutions. 

As fa/N increases even further, we move into the range where Tg < 10 MeV, and 
saxion decay might disrupt BBN. In this excluded region, two solutions to the restriction 
^adh'^ = 0.1123 appear. The first has dark matter dominated by thermal axinos and ultra- 
high Tr S> fa/N wherein the axinos and axions are severely diluted by saxion entropy 
production; these solutions are not exhibited on the plot. The second solution allows for 
much lower Tr values in which case dark matter is dominated by axion production, albeit 
with some dilution due to coherent oscillation production of saxions. These high fa/N 
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Figure 12: Upper frame a): Contribution of axions and TP and NTP axinos to tire DM density as 
a function of the PQ breaking scale fa/N, for ^^^h? = 10 and m^^ — 430 GeV, with ma — 100 keV 
and 6i — 0.05; is adjusted such that ilaah^ — 0.1123. We assume = ttlq ~ 1 TeV and 
tjIq = Lower frame b): the value oiTu (and Tr/t) that is needed to achieve ilaah^ = 0.1123 

for nia — 0-1 and 1 MeV. 

solutions, however intriguing, are all excluded because such high values of /a/A^ suppress 
the saxion and Zi lifetimes, so their decays will affect BBN. 

The upshot of Fig. || is that, for fa/N slightly below 10^^ GeV, the value of Tr 
has increased to over 10^^ GeV while maintaining Qaa.h'^ = 0.1123. Although the saxion 
entropy injection leads to higher values of Tr (when compared to Fig.|l|), the allowed range 
for the relevant temperature for leptogenesis (T^/r) is actually reduced, due to the BBN 
bounds on the decaying saxion. For the case of rria = 0.1 MeV, the value of Tff/r reaches 
to ~ 10^'^ GeV for fa/N ~ 10^'^ GeV, thus potentially reconciling thermal leptogenesis 
with the gravitino problem in the AY scenario. However, the solution with = 1 MeV 
never reaches as high as Tfj/r ~ 2 x 10^ GeV before entering the BBN-excluded region, and 



-19- 




so in this case does not lead to a reconciliation of themal leptogenesis with the gravitino 
problem. Comparing Figs. || and |l^ we can see that, at least for this case, the range of fa 
values which accomodates thermal leptogenesis is actually reduced when the saxion entropy 
injection effect is included. 

3.3 Scan over PQMSSM parameters including dilution due to entropy injection 
from saxion decay 



While Fig. 12 holds for particular values of the PQMSSM parameters, save for fa/N , we will 



now scan over the remaining PQ parameters rria and Oi, as well as fa/N, as in Section 'Lb 
This time, we will adopt = mo = 1 TeV and Si/{fa/N) = 1, and allow for saxion- 
induced entropy dilution of mixed axion/axino DM according to the procedure described 



in the last section. The results for h'^ = 0.04 and = 430 GeV are shown in Fig. 13 
where we plot the value of Tfj/r needed to maintain ^aah"^ = 0.1123 versus PQ breaking 
scale fa/N. The line where r = Sf / St = 1 is shown in magenta. The red points violate 
BBN bounds due to late decaying Zi, while the green points violate BBN bounds due to 
late-time saxion decays (T^ < 10 MeV). The light blue points have > 20% WDM or > 1% 
HDM, while the dark blue points satisfy all constraints. The WDM/CDM bound following 
Boyarsky [37] is again indicated as a dashed blue line. By including dilution of DM from 
saxion production and decay, the allowed points can reach to Tr as high as ~ 10^'^ GeV for 
fa/N ~ 1.5 X 10^'^ GeV although the value of Tfj/r reaches only as high as ~ 10^^ GeV. 
These points with Tji/r ^ 2 x 10^ GeV evidently reconcile thermal leptogenesis with the 
gravitino problem even in the presence of entropy injection from saxion decay. 

In Fig. |l^, we plot the axion mis- alignment angle 9i. Unlike the previous results in 
Fig. 1^ with no entropy injection, the allowed values of 9i with T^/r > 2 x 10^ GeV span a 
range from to ~ 1.5 radians: for higher values of Tr, larger values of 9i can be tolerated 
since the relic abundance of axions is now diluted by saxion decay. Therefore, in this case. 
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Figure 14: Allowed and disallowed points in the Oi vs. Tr/v plane for fl^ = 0-04 and 
430 GeV, with m, = 1 TeV. 



the axion mis-alignment angle is not required to take unnaturally small values, as opposed 
to the case without the saxion dilution, shown in Fig.^. 

To see whether axinos or axions dominate the DM density including entropy from 
saxions, in Fig. ^ we plot the same points, but this time versus axion relic density ila/i^. 
We see that the bulk of points with T/j/r > 2 x 10^ GeV that are BBN-allowed indeed have 



mainly axion CDM. Note that the point shown in Fig. 12, which has Oi = 0.05 and mainly 
axino DM (at Tr > 2 x 10^ GeV), corresponds to the few points of Fig. 15 at low ri^/i^ and 
is not the most common scenario, since it requires quite small values of the mis-alignment 
angle. Given that fa/^ ~ 3 — 15 x 10^^ GeV, we expect the axion mass rua ~ 0.4 — 2 //eV, 
somewhat below the range where AD MX is searching [49]. 

3.4 More general scan over MSSM parameters 

Next, we generalize our results for a general PQMSSM model, where we now allow and 
to be free parameters included in our scan, as in Fig. For simplicity we keep the 
saxion mass fixed at nis = 1 TeV. We keep only points with T^/r > 2 x 10^ GeV, which 
potentially reconcile thermal leptogenesis with the gravitino problem. The result is shown 
in Fig. ^ where the red points are excluded due to the BBN constraints on Zi decays; no 
green points due to constraints from BBN on saxion decay are visible. By comparing Figs. 
|6| and 16 , we see that due to the saxion dilution of the neutralino relic density, the BBN 



bounds on are less severe and a larger portion of the MSSM parameter space can be 
consistent with thermal leptogenesis. 

3.5 Thermal leptogenesis-allowed regions of the mSUGRA plane 

As a last point of this study, let us apply our general results to the showcase mSUGRA 
model in the mo vs. mi/2 plane. In order to make our results independent of a particular 
choice of PQ parameters, we consider the bounds on ^^i^^ obtained from the general 
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Figure 16: Allowed and disallowed points in the ^^i^^ vs.ra^^ plane for a general scan over SUSY 
models with a bino Zi and nis — 1 TeV. For all points, we require T^/r > 2 x 10^ GeV. In this 
plot, we include the effect of entropy production from saxion decay. Dark blue points are consistent 
with BBN and have mainly CDM with at most 20% WDM and/or 1% HDM admixture. The region 
below the dashed line represents the MSSM parameter space where 99% of the DM/BBN consistent 
solutions lie when applying weaker WDM/HDM requirements as discussed in the text. 



PQMSSM scan in Sections 2^ and |3.4| (for the case with saxion entropy injection). These 
bounds are represented by the dashed hnes in Figs. |6| and We may then translate this 
into a contour in the mo vs. mi/2 plane of mSUGRA for y4o = 0, /U > and constant 
tan/3, as shown for the cases of tan/3 = 10, 50 and 55 in Fig. 17. The gray regions are 
excluded because they violate the LEP2 limits on Higgs or sparticle masses^ or have a 



^The LEP2 limit on a SM-like Higgs scalar h is > 114.4 GeV. Here, we use mu > 111 GeV allowing 
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stau as next-to-next-to-lightest SUSY particle (NNLSP); this latter case requires special 
treatment as for example in Ref. [51]. 

In frame a), we show the mSUGRA mo vs. mi/2 plane for tan/3 = 10. The strips 
of dark blue and purple points show the regions that allow for Tr > T™" = 2 x 10^ 
GeV, while maintaining O.aa.h'^ = 0.1123 and respecting bounds from BBN. The subset of 
purple points at low mi/2 satisfies in addition the following constraints on low energy (LE) 
observables: 

1. Aa^^^^ = (7.90 - 37.39) x 10"^° , 

2. BR{b S7) = (2.79 - 4.3) x 10-^ 

3. BR{Bs fi+i-r) < 4.7 X 10-^ 

4. 0.55 < BR{Bu T+Ur)^^^^/BR{Bu ^ r+i^rf^ < 2.71 

where 1. — 3. were calculated using Isajet/Isatools and 4. was calculated using Superlso. 

We see that the AY consistent regions, although broader, are very similiar to the classic 
mSUGRA regions with neutralino dark matter: the stau co-annihilation region at low mo 
and the light Higgs resonance region where ZiZi —^h at mi/2 ~ 150 GeV. The reason is 
that a rather low abundance of thermal neutralinos is required in the AY scenario to satisfy 
BBN constraints on late decaying Zis. For comparison, the classic mSUGRA strips where 
the neutralino relic density ^^_^h'^ = 0.1123 ± 0.0105 are shown as yellow/orange points. 

Invoking next the vs. m^^ contour of Fig. which includes the effect of entropy 
generation from a rus = 1 TeV saxion, the AY-consistent regions broaden out considerably. 
The region with Tfj/r > T™" is denoted here by light blue points, and expands to fill the 
lower mo portion of the mo vs. mi^2 plane along with a band around mi^2 ~ 400 where 
turn-on of the ZiZi — >■ tt annihilation chanel reduces the neutralino abundance. The 
portion of the leptogenesis consistent region including saxion decays and LE constraints 
is colored in pink, and requires mi/2 < 550 GeV and mo < 500 GeV, so as to allow for 
a significant contribution to {g — 2)^ by light charginos and sneutrinos. The remaining 
unshaded (white) region of the mSUGRA plane does not allow for an AY reconciliation 
of thermal leptogenesis with the gravitino problem, with or without saxion decays, mainly 
because the relic density of neutralinos is so large that the BBN constraints on late decaying 
Zi are violated. 



Frame b) of Fig. 17, shows the analogous plot for tan/3 = 50. In this case, b- and 
r- Yukawa couplings increase greatly, while the value of m^ drops, enabling efficient an- 
nihilation of neutralinos via stau coannihilation or s-channel A exchange. The apparent 
neutralino abundance ^^_^h'^ is severely reduced, and less constrained by BBN. The area of 
leptogenesis-consistent regions increases. Furthermore, the SUSY contributions to 6 ^ S7 
and {g — 2)^ increase with increasing tan/3, and so the region which is consistent with LE 
constraints moves to higher mi/2 values. If saxion entropy production is added, almost the 
whole plane is allowed by the AY scenario. 



for an approximate 3 GeV error on the theory calculation of rrih. For the SUSY mass limits we use those 
implemented in Superlso [50]. 
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Figure 17: Regions in the mo vs. miji plane of the mSUGRA model with = and > 
which satisfy 1. Tr > T^™ = 2 x 10^ GeV (dark blue), 2. Tr > T™™ and LE constraints (purple), 
3. Tr/t > r™™ with saxion entropy injection (light blue) and 4. Tr/t > TJ^™ with saxion entropy 
injection and LE constraints (pink). For comparison, the yellow/orange points indicate the classic 
mSUGRA regions with ft^^h'^ = 0.1123 ± 0.0105. We show frames for a) tan/3 = 10, b) tan/3 50 
and c) tan /3 — 55. 

Finally, frame c) shows the case of tan /3 = 55, where the ^-resonance dominates 
the ZiZi annihilation amplitudes. Here, we see that a huge swath of parameter space 
is AY-consistent, even without the effect of saxion decays. By including entropy from 
saxion decay, the entire mo vs. mi/2 plane becomes AY-consistent. The part which is 
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consistent with LE constraints follows suit, leading to a large region of parameter space 
that is consistent with all constraints. 

4. Conclusions 

In this paper, we reported on investigations of the viability of the Asaka-Yanagida sug- 
gestion that a mass hierarchy with m(sparticle) > m-g > can be used to reconcile 
thermal leptogenesis, which requires Tr > 2 x 10^ GeV, with the gravitino problem, 
which seemingly requires much lower Tr to avoid BBN constraints and overproduction 
of neutralino dark matter. In the AY scenario, the G decays inertly to aa. BBN con- 
straints on sparticle ^ G -|- particle can be avoided because the much faster decays 
sparticle a + particle are now allowed. We re-examined the AY scenario in Sec. § by 
including 1. updated measurements on the total dark matter abundance Q.^Aih'^ — 0.1123, 
2. updated calculations of thermal axino and gravitino production, 3. the contribution of 
relic axions and 4. BBN constraints on late decaying Zis. Furthermore, in Sec. |3|, we in- 
cluded dilution of dark matter by saxion production and decay. The latter effect can be 
neglected if rris is in the multi-TeV range and the initial saxion field strength Sj is somewhat 
smaller than the PQ breaking scale fa/N. 

We found in Sec. |2|, neglecting the saxion entropy effect, that the AY scenario does 
work under the conditions that (i) fa/^ is rather large > 10^^ GeV, implying a somewhat 
lighter axion than is presently searched for by ADMX [49], (ii) the apparent neutralino 
relic density ^^_^h'^ is not too big: ^^_^h'^ ^ 1, (Hi) the value of is at least in the several 
hundred GeV range in order to hasten the Zi decay rate, and (iv) the axion mis-alignment 
angle 6i is on the small side < 0.5 to suppress overproduction of axions when fa/N is large. 

By including saxion production and decay in Sec. ^, we can dilute the axino and 
also axion DM abundance, which in turn allows for somewhat higher values of Tr up to 
~ 10^'^ GeV to be generated. However, since saxion decay also dilutes the baryon density, 
in this case we must require instead Tji/r > 2 x 10^ GeV. The saxion mass rrig needs to 
be rather large to avoid BBN constraints on late decaying saxions if Tr is to be high. In 
this case, the DM is likely to be mainly axions, although a few cases with mainly axino 
DM were generated. The axion mis-alignment angle need not be small here since the axion 
abundance can be suppressed by entropy injection from saxions. We have also found that 
a large portion of the MSSM parameter space {^^_^ and rn^^) can be consistent with high 
Tr and still avoid the BBN bounds on late decaying neutralinos, due to the dilution of the 
neutralino relic density after the entropy injection from saxion decays. 

The observable consequences of our final results are as follows. If the AY scenario with 
m(sparticle) > > is to reconcile thermal leptogenesis with the gravitino problem, 
then we expect several broad results to ensue: 

1. discovery of SUSY at the LHC, with a reconstructed ^^i^'^ large, lest Zis 

are produced at too large a rate in the early universe, and their late decays disrupt 
BBN; 
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2. a SUSY mass spectrum consistent with SUGRA models with a rather hght (but still 
weak scale) gravitino, since the gravitino mass must be lighter than all observable 
sparticles; 

3. a mainly bino-like Zi, to quicken decays into a^y/Z^ with mass m^^ in the hundreds 
of GeV range, which also helps diminish the lifetime; 

4. no direct or indirect detection of neutralino (WIMP) dark matter; 

5. finally, we expect discovery of an axion to be likely, but in the mass range ~ 0.1 — 2 
//eV, somewhat below the values presently being explored. 
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